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La4Cu3MoOi2 is a cluster antiferromagnet where copper spin-1/2 form a network of strongly 
coupled spin-trimers. The magnetic properties of this material have been examined using magnetic 
neutron scattering. At low temperatures, excitations from the ground state are observed at 7.5(3) 
meV and 132.5(5) meV. An additional peak in the neutron scattering spectrum, which appears at 
125.0(5) meV on heating is ascribed to a transition between excited states. The wave-vector and 
temperature-dependence of the inelastic magnetic scattering cross section is consistent with intra- 
trimer transitions. Magnetic neutron diffraction reveals antiferromagnetic order below Tat=2.6 
K that doubles the unit cell along the a direction. The ordered magnetic structure is described 
as inter-trimer order where spin correlations within trimers are controlled by the strong intra- 
trimer interactions. Combining the information derived from elastic and inelastic magnetic neutron 
scattering with group theoretical analysis, a consistent set of intra-trimer interactions and ordered 
magnetic structures are derived. The experiment provides a simple worked example of magnetism 
associated with inter-atomic composite degrees of freedom in the extreme quantum limit. 

PACS numbers: 75.10.Jm, 75.25.+Z 



I. INTRODUCTION 

Geometrically frustrated magnets are distinguished by 
an anomalous cooperative paramagnetic phase extending 
to temperatures well below characteristic microscopic en- 
ergy scalesE Fluctuations in this phase are strong and 
non-trivial in that they satisfy certain local constraints. 
In some cases, further cooling fails to produce a phase 
transition _aricl the low temperature state is quantum 
disorderecflu. In other cases, however, thermal or quan- 
tum fluctuations, magneto-elastic coupling, impurities, 
or sub-leading exchange interactions yjftlrirp finite value 
of (S) at sufficiently low temperaturesffiljLl. 

To advance the understanding of geometrically frus- 
trated magnetism, we have explored spin correlations 
in La4Cu3MoOi2, a material which affords a spectacu- 
lar example of geometrical frustration in the quantum 
limit. La4Cu3MoOi2 is a monoclinic ABO3 type cupratej 
with space group P112i/m and B-cations in 1:3 ratio.tll 
Shown in Fig. [l], the CU3M0O4 layers of the material 
are built from CU3O triangular clusters, which we shall 
denote trimers. Consideration of the crystal structure 
and the Goodenough rulesO indicate that intra-trimer 
exchange interactions are orders of magnitude stronger 



than-jfttcr-trimer interactions. Magnetic susceptibility 
dataEJ'O show two distinct linear regimes for l/x(T) ver- 
sus temperature. For T > 400 K the Curie Weiss fit 
yields ^ e //=1.81 \ib and &w=— 558 K, which is consis- 
tent with strong intra-trimer interactions. For T < 250 
K the Curie constant decreases by a factor of 0.39 and 
the Weiss constant decreases to —16 K. This indicates a 
cross over to a cooperative paramagnetic phase where 
each trimer represents a composite spin-1/2 degree of 
freedom. The small Weiss temperature in this phase in- 
dicates weak AFM inter-trimer interactions. At lower 
temperatures still, a maximum in x(^) an d a peak in the 
specific heat indicates_an antiferromagnetic phase tran- 
sition at T N = 2.6 K.EJ 

In this paper, we present a comprehensive neutron 
scattering study of La4Cu3MoOi2- Polarized and unpo- 
larized inelastic neutron scattering data provide evidence 
for magnetic excited states 7.5 meV and 132.5 meV above 
the ground state. A calculation based on a simple intra- 
trimer Hamiltonian is presented to explain the energy 
levels, the wave vector dependence, and the temperature 
dependence of the inelastic neutron scattering cross sec- 
tion. Using elastic neutron diffraction we also show that 
La4Cu3MoOi2 has long range antiferromagnetic order at 
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FIG. 1: CU3M0O4 plane of La4Cu3MoOi2. Actual coordi- 
nates for copper atoms are given in appendix A. The solid 
lines represent the triangle clusters. Ji, J2 and J3 are the in- 
tratriangle couplings; J's are the weak intertriangle couplings 
that are assumed to be antiferromagnetic and of similar ma* 
nitude as those in the mean field analysis of Wessel and Haas.u 
The arrows illustrate one plane of the T2 or T3 magnetic struc- 
tures with -i/i = listed in table O. 



low temperatures. While a unique ordered structure can- 
not be identified directly from the diffraction data, anal- 
ysis based on ordering of composite spin- 1/2 degrees of 
freedom on neighboring trimers yields a consistent set of 
intra-trimer exchange interactions and inter-trimer spin 
configurations. The experiments and analysis provide a 
comprehensive understanding of geometrically frustrated 
quantum magnetism in a simple model system with im- 
portant analogies to more complex systems. 



II. EXPERIMENTAL RESULTS 

A powder sample of La^uaMoCUa was synthesized us- 
ing a previously published method £j Rietveld analysis of 
powder neutron diffraction data confirmed a single phase 
sample with space group P112i/m and lattice parame- 
ters a = 7.9119(5) A, b = 6.8588(4) A, c = 10.9713(3) A, 
and 7 = 90.008(9)° at 10K. The elastic neutron scatter- 
ing measurements were carried out on the BT2 thermal 
neutron triple-axis spectrometer at the NIST Center for 
Neutron Research. For that experiment we used a 30.2 g 
powder sample in a cylindrical container with a diameter 
of 1.6 cm. Pyrolytic Graphite (PG) crystals set for the 
(002) reflection were used to select 14.7 meV neutrons for 
diffraction. There was a PG filter in the incident beam 



to suppress higher order contamination and collimations 
were 60' — 40' — 40' — 200' through the instrument from 
source to detector. 

Inelastic neutron scattering measurements were per- 
formed on the HET direct geometry time-of-flight spec- 
trometer at ISIS pulsed spallation neutron source. For 
that experiment we used a 80.2 g powder sample in a 
cylindrical container with a diameter of 3.1 cm. We used 
that instrument's "sloppy" chopper at frequencies of rev- 
olution 250 Hz and 400 Hz for incident energies of 40 meV 
and 160 meV respectively. 

Polarized inelastic neutron scattering measurements 
were performed on the C5 triple-axis spectrometer at 
the NRU reactor in Chalk River Laboratories in Canada. 
Magnetized Heusler crystals set for the polarizing (111) 
reflection were used as monochromator and analyzer with 
the latter fixed to reflect 14.6 meV neutrons. A cold 
sapphire filter was placed before the monochromator to 
eliminate high-energy neutrons. A PG filter was placed 
in the scattered beam to suppress order contamination 
at the analyzer. An energy dependent correction was ap- 
plied to the incident monitor count rate to account for 
incident beam A/2 contamination. For this experiment 
we used a 120 g powder sample in a cylindrical container 
with a diameter of 2.5 cm. The collimation was 45' in 
the incident and 80' in the scattered beam. A Mezci flip- 
per was placed in the scattered beam. The flipping ratio 
measured at 14.6 meV on a powder reflection was 24:1. 

Absolute normalization turned out to be crucial to 
derive information about the magnetic structure. We 
used nuclear Bragg scattering to normalize the magnetic 
Bragg peak intensity!] The inelastic neutron scattering 
data was normalized by comparison to incoherent elas- 
tic count rates from a 20.14 g vanadium sample in the 
ISIS experiment. A correction for neutron absorption in 
La4Cu3MoOi2 was applied in this latter case. 



A. Elastic Neutron Scattering 

Fig. ^| shows the temperature dependence of mag- 
netic neutron diffraction at Q=0.397 A^ 1 . The data 
is evidence for an antiferromagnetic phase transition at 
Tat=2.6 K <C Qw- The inset shows temperature dif- 
ference data indicating that elastic scattering from the 
low temperature spin structure is in the form of a well- 
defined Bragg peak. To within error the location of the 
peak is at ^a* implying that the magnetic order doubles 
the unit cell in the a direction. After due consideration 
of the finite instrumental resolution (solid bar in the in- 
set), we derive a lower limit of 100 A for the magnetic 
correlation length. 

B. Inelastic Neutron Scattering 

Fig. H shows the energy dependence of the scattering 
cross section at various temperatures. Because of the 
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FIG. 2: Temperature dependence of the (|00) magnetic 
Bragg peak intensity. The inset shows the wave vector de- 
pendence of the difference between elastic magnetic neutron 
scattering at T = 0.4 K and T = 10 K. The energy resolution 
was 1.3 meV and the wave vector resolution is indicated by 
the solid line in the inset. 



wide dynamic range, we show only the interesting mag- 
netic parts of the spectrum, low energies on the left and 
higher energies on the right. The right panels clearly 
show an excitation at 132.5 meV. As temperature in- 
creases, the peak intensity at 132.5 meV decreases and 
a second peak emerges at Ticu—125 meV. A likely expla- 
nation for the lower energy peak is that it corresponds 
to a transition to the 132.5 meV state from an excited 
state at 132.5 meV-125 meV=7.5 meV, which is popu- 
lated on heating. To explore this scenario, the left panels 
focus on the energy range around 7.5 meV. As is ap- 
parent from the increase of intensity with temperature, 
there is significant phonon scattering in the lower energy 
range and this complicates the task of isolating the mag- 
netic contribution to the scattering cross section. Since 
phonon scattering increases with temperaturec3 while low 
energy magnetic scattering generally decreases with tem- 
perature, we used high temperature data (T = 200 K) 
where magnetic scattering is negligible to determine the 
phonon density of states. At each temperature, the ap- 
propriate thermal factor was then applied to yield the 
phonon contribution to inelastic neutron scattering. The 
triangles in the left panels of Fig. || show the background 
subtracted data that indicate a residual peak in the ex- 
citation spectrum, which we tentatively associate with 
magnetic scattering. 

To unambiguously verify the existence of a magnetic 
excitation at 7.5 meV, a polarized neutron scattering 
measurement was carried out at T=6 K and at a mo- 
mentum transfer of 1.5 A" 1 . Only magnetic scattering 
can produce a difference between spin flip scattering with 
guide fields at the sample position parallel and perpen- 
dicular to wave vector transfer Ii3 Such difference data 
is shown as closed triangles in Fig. |^(a) and they pro- 
vide immutable evidence for a magnetic excitation at 8(1) 
meV. Also shown is the sum of spin flip and non-spin-flip 
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FIG. 3: ftoj-dependence of the normalized Q-averaged neu- 
tron scattering intensity I(Q,hu>) at various temeratures and 
in two different energy ranges. The open circles are raw ex- 
perimental data. The open triangles are the phonon back- 
ground subtracted data in (a)-(c), and the excess scattering 
above that at 10 K for (e)-(f). The closed circles and closed 
triangles in (a) are total scattering and magnetic scattering 
intensities respectively from the polarized neutron measure- 
ment. A single scale factor was applied to the polarized data 
for best agreement with the unpolarized time of flight data. 



scattering as solid circles. Twice the count rate with the 
analyzer turned through 4 degrees was subtracted from 
the summed data and a single overall scale factor was 
applied to all the polarized data to facilitate comparison 
to the time of flight data (open circles). The excellent 
agreement between the two independent determinations 
of the magnetic contribution to inelastic scattering at 7.5 
meV provides strong evidence for an intra-trimer excited 
state at this energy. 

From gaussian fits to all inelastic data, we derived 
the temperature dependence of the integrated intensi- 
ties, peak positions and peak widths that are reported 
in Fig. [|. For peak widths, the calculated energy res- 
olution was subtracted in quadrature to produce values 
for the intrinsic half width at half maximum relaxation 
rate, T or inter-trimer bandwidth. The energy levels are 
temperature independent to within error as expected for 
intra-trimer excitations. T is of order the low-T Curie- 
Weiss temperature, which is consistent with inter-trimer 
interactions being the main source of intra-trimer level 
broadening. Transitions involving the excited state dou- 
blet appear to have a greater relaxation rate than those 
involving other levels. Possible reasons for this include 
magneto-elastic effects and enhanced inter-trimer cou- 
pling for trimers occupying the 7.5 meV excited state. 
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FIG. 4: Temperature dependence of parameters char- 
acterizing high and low energy magnetic excitations in 
La4Cu3MoOi2. Frame (a) shows the integrated intensities 
for the 132.5 meV (open circles), the 125 meV (triangles), 
and the 7.5 meV (closed circles) modes. Frame (b) shows the 
peak positions and frame (c) shows the intrinsic half width at 
half maximum relaxation rate for each mode. Lines in frame 
(a) were calculated from the trimer model. 
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FIG. 5: Wave- vector dependence of energy integrated inten- 
sity for the 132.5 meV (circles) and the 7.5 meV excitations 
(triangles) in La4Cu3MoOi2. Solid lines show the calculated 
Q-dependence of the neutron scattering cross section from Eq. 
(p36|) and (B7) with no adjustable parameters. 



trimer have the symmetry of an equilateral triangle there 
is a fourfold degenerate ground state composed of two 
degenerate Kramers doublets. For lower symmetry the 
degeneracy is lifted and there are two low energy dou- 
blets. While they approximate the symmetry of isosceles 
triangles, the spin- 1/2 trimers in La4Cu3MoOi2 have no 
exact symmetry elements. In Appendix B, we calculate 
the energy levels and scattering cross sections for a gen- 
eral spin-1/2 trimer with a model Hamiltonian of the 
form 



H — JiSi • S2 + J2S2 ■ S3 4- J3S3 • Si 



(1) 



Fig. g shows the wave- vector dependence of the energy 
integrated intensity for the 7.5 meV and the 132.5 meV 
modes. Kinematical limitations prevented measurement 
of the high energy excitations over a significant range of 
wave-vector transfer and at lower energies admixture of 
phonon scattering complicates the analysis. Still there is 
evidence that the intensity of the 7.5 meV mode decreases 
in rough correspondence with the magnetic form factor 
for copper, which indicates that this excitation involves 
a very small cluster of copper atoms. In addition we 
shall see that the relative intensity of the two modes is 
consistent with a simple trimer exchange Hamiltonian. 



III. ANALYSIS AND DISCUSSION 

A. Energy Level Scheme 

From the energy level scheme we can extract de- 
tailed information about intra-trimer exchange interac- 
tions. The highest energy state for a spin-1/2 trimer 
with antiferromagnetic Heisenberg interactions is a quar- 
tet with total spin 3/2. If the interactions within the 



The splitting between the d oublets is 
Eoi = v^i 4- Jf 4- Jf - Ji J2 - J2J3 ~ J3J1 and the 
splitting between the ground state doublet and the quar- 
tet is E 2=^(Ji 4- J 2 4- J 3 + Eqi). These energies can be 
associated respectively with the 7.5 meV and the 132.5 
meV transitions observed by inelastic neutron scattering. 
From this we can derive the average intra-trimer ex- 
change constant to be J = {2Eq2 — Eqi)/ 3 = 85.8 meV. 
The splitting between doublets yields information about 
the ratios between intra-trimer exchange constants. Ra- 
tios that are consistent with the 7.5 meV doublet-doublet 
transition lie on an ellipse in the J1/J3 versus J2/J3 plot 
shown in Fig. ^ The ellipse is centered at the equilateral 
point (1,1) and the major axis lies along the isosceles 
J 1 = J 2 line. The half major axis is approximately 
V / 9E% 1 /(2E% 2 - 2E al E Q2 ) = 0.124 and the half minor 
axis is approximately ^3E^ 1 /(2E^ 2 - 2E Q iEq 2 ) — 0.071. 
In the following, we shall label intra-trimer interaction 
parameters that are consistent with the spectroscopic 
information by the counter-clockwise azimuthal angle, 
ip, on this ellipse with tp = corresponding to J\ = J 2 
and ip = mr describing isosceles triangles. 

Wessel and HassB recently studied the phase diagram 
for La4Cu3MoOi2 using a model Hamiltonian including 



5 



2.0 



1.5 



1.0 



0.5 - 



0.0 



1 


(io) 


(o,i) 


(ii 


\ 

1 


(io) 

(o.-l) 

1 



0.0 



0.5 



1.0 

JiA 



1.5 



2.0 



FIG. 6: Zero-temperature magnetic phase diagram of the 
AFM spin- 1/2 trimer square lattice with a weak intertrian- 
gle coupling J'=0.01J3. The excitation energies observed in 
the present experiment imply that the ratios of exchange con- 
stants lie on the ellipse shown in the center of the figure. The 
observation of (^00) magnetic order also helps to constrain 
possible values of the exchange constants. 



the nearest neighbor interactions within the a-b plane 
defined in Fig. H. This study shows that different ratios 
of intra-trimer interactions yield different in-plane wave 
vectors for long range order. The wave vectors predicted 
in the limit of vanishing intcr-trimer interactions at T = 
are indicated on Fig. [| Given that Tjy <C J this is 
the appropriate limit to consider. Taking into account 
this T = phase diagram leads to the conclusion that 
the azimuthal angle specifying intra-trimer interaction 
asymmetry must satisfy ip = 7r or — 7r/4 < ip < tt/4. 

In Fig. I the wave vector dependence of the energy 
integrated intensities at T — 10 K are compared to the 
formulae of Eq. (B6) and Eq. (B7). The agreement be- 
tween model and data is quite satisfactory considering 
that there are no adjustable parameters. While the com- 
parison does not provide information on intra-trimer ex- 
change, it supports the identification of a magnetic con- 
tribution to inelastic scattering at 7.5 meV, as the abso- 
lute cross section inferred from polarized and unpolarizcd 
magnetic neutron scattering is in perfect agreement with 
that predicted for the transition between the two dou- 
blets of the spin-trimer. 

The temperature dependence of the integrated intensi- 
ties for t he three inelastic peaks follows from Eq. (B6) to 
Eq. (B8). It depends only on the population of the intra- 
trimer levels, which in turn only depends on Eq\ and 
E02 ■ The solid and dashed lines in Fig. [| (a) were calcu- 
lated from these formulae using the experimental values 
for E i and Eq 2 and they are found to be in excellent 
agreement with the data. 



TABLE I: The basis functions for irreducible representations 
of space group P112i/m with magnetic wave vector (iOO) and 
for atoms on the 2e site. The exact atomic coordinates are 
given in appendix A. 
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B. Magnetic Structure 

Only the (5 00) magnetic Bragg peak could be detected 
in the present experiment. We obtained the absolute 
magnetic structure factor at (^00) by comparison to the 
(230) nuclear Bragg peak. By combining this informa- 
tion with the symmetry analysis described in Appendix 
A, and the spectroscopic information presented in sec- 
tion [III A , we can associate each of a few possible or- 
dered structures with a specific set of intra-trimer ex- 
change constants. 

The symmetry analysis is based on the assumption 
that any magnetic structure adopted through a second 
order phase transition can be expanded in basis func- 
tions for a single irreducible representation of the mag- 
netic space group. Table | lists the basis functions for four 
of the irreducible representations of magnetically ordered 
La4Cu3MoOi2- ti and T4 describe uniaxial spin configu- 
rations with spins oriented along the c direction while t 2 
and T3 correspond to co-planar structures with spins in 
the a — b plane. 

To make the connection between spectroscopy and 
magnetic structure we consider the ordered structure not 
as ordering of individual spins, but as ordering of the 
composite spin-1/2 degree of freedom associated with the 
lowest energy intra-trimer doublet. Appendix B lists the 
eigenstates of the intra-trimer spin Hamiltonian for ar- 
bitrary exchange interactions. Neel order corresponds to 




FIG. 7: (a) Expectation value for the magnetic moment along 
a small applied field and at T — for the three copper ions in a 
trimer as a function of index angle ip. ip is the azimuthal angle 
spanning intra-trimer exchange interactions that are consis- 
tent with a 7.5 meV doublet-doublet transition (see Fig. |^). 
J1/J3 is equal to J2/J3 at ip = and n. (b) Elastic magnetic 
neutron scattering cross section for the (^00) Bragg peak per 
crystal unit cell calculated as a function of ip when the spins 
are along c direction and adopt the magnetic structure as- 
sociated with the T4 irreducible representation. Dashed line 
shows the actual measured magnetic Bragg intensity. 



choosing a specific quantization axis on each spin trimer 
and alternating the doublet occupation consistent with 
the (^00) magnetic wave- vector. Fig. ^|(a) shows the spin 
projection on the quantization axis, for each of the three 
atoms on a trimer as a function of the azimuthal angle, 
ip, that indexes possible intra-trimer exchange constants. 
Given a quantization axis, and a value for ip, the mag- 
netic structure factor for the (^00) Bragg peak can be 
calculated for comparison with the measured absolute 
intensity of the magnetic Bragg peak. 

Spin structures corresponding to irreducible represen- 
tations 7~i are inconsistent with the measured magnetic 
Bragg peak intensity for all values of ip. For T4 Fig. (?](b) 
shows that the calculated intensity is consistent with 
the measured intensity, for four different values of ip — 
(0.01(2), 0.32(4), 0.64(4), 0.94(2))tt. Of these solutions, 
only the former two are consistent with the stability anal- 
ysis for the (|00) structured and of these only the first 
is consistent with the approximately isosceles spin trian- 
gles. 

Irreducible representations Ti and T3 describe spin con- 
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FIG. 8: Combination of index angle ip and the spin orientation 
angle <p for magnetic structures with irreducible representa- 
tion T2 or T3 that yields a (^ 00) magnetic Bragg intensity 
consistent with the experiment. Parameters in the hatched 
areas stabilize theJ^OO) magnetic structure for weak inter- 
trimer interactionsa. 



TABLE II: Possible spin configurations corresponding to 
isosceles triangles with Ji = J2. The exact atomic coordi- 
nates are given in appendix A. 
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figurations with spins in the a-b plane. For simplicity we 
consider only uniaxial spin configurations spanned by the 
angle <p between the spin-trimer quantization axis and the 
a direction. T2 and T3 are different only in the registry of 
magnetic order in different a — b planes. It can be shown 
that T2 and T3 cannot be distinguished so they must be 
considered simultaneously. Fig. || shows values for ip and 
<f> that are consistent with the measured (^00) magnetic 
Bragg intensity. Also indicated as hatched areas are the 
values of ip where the (|00) structure is stable.cl 

As previously mentioned the trimers have an approx- 
imate mirror plane, and this leads to the expectation 
that Ji w J2. The seven spin structures and values of 
exchange constants that are consistent with this assump- 
tions are listed in table ||. Of these, the three structures 
with ip = seem more likely because as opposed to the 
ip = 7T structures, they have a finite range of stability in 
the J1/J3 — J2/J3 plane (see Fig. ||). 
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IV. CONCLUSION 



APPENDIX A: SYMMETRY ANALYSIS OF 
MAGNETIC ORDERING 



In summary we have presented neutron scattering data 
that provide detailed microscopic information about frus- 
trated quantum magnetism in La4Cu3MoOi2. Neu- 
tron spectroscopy yields the average intra-trimer ex- 
change constant of 85.8 meV and narrows possible intra- 
trimcr exchange ratios to an elliptical trajectory in the 
J1/J3 ~ J2/J3 plane. Neutron diffraction provides ev- 
idence for low temperature inter-trimer magnetic order 
that doubles the unit cell along the a direction. The 
(^00) magnetic Bragg peak intensity and the approxi- 
mate isosceles nature of the spin triangles narrows the 
possible exchange constants to J% = J2 = 88.2 meV and 
J 3 = 81.1 meV or J 1 = J 2 = 83.2 meV and J 3 = 91.1 
meV. A previously published mean field analysis of mag- 
netic order in La4Cu3MoOi2 indicates that the former 
combination of exchange constants is most likely. The 
corresponding ordered spin structures are uniaxial with 
two parallel spins at the base of the isosceles triangle 
((S z ) = 3-) and an antiferromagnetically correlated spin 
of half the magnitude at the apex. The data are consis- 
tent with a spin direction either along b or c. Another 
unresolved issue is the magnetic stacking sequence along 
c for which there are two options for spins oriented along 
b. 

There is an instructive analogy between the frus- 
trated cluster antiferromagnetism in La4Cu3MoOi2 and 
rare earth magnets. In materials with Kramers rare 
earth ions, intra-atomic correlations establish effective 
spin-1/2 degrees of freedom, which subsequently develop 
long range magnetic order due to inter-atomic exchange 
interactions .til Magnetic neutron scattering from_such a 
system carries the rare earth atomic form factor.E3 Mag- 
netism in La 4 Cu 3 MoOi2 is also based on a composite 
spin-1/2 degree of freedom, only it is spread over three 
atoms, it is established by inter-atomic exchange interac- 
tions, and it carries an oscillatory trimer "form factor" . 
La4Cu3MoOi2 is thus a particularly simple example of a 
concept of increasing importance in quantum piag netism. 
End states and holes in Haldane spin chains ,li§li3 r impu- 
rity spins in high temperature superconductors ,113 and 
spontaneously formed or structurally defined spin clus- 
ters in frustrated magnets. t§Ej All are strongly corre- 
lated systems where suitably defined multi-atom com- 
posite spin degrees of freedom provide an enormous sim- 
plification for understanding low energy spin dynamics 
and the corresponding thermodynamic properties. 



A magnetic structure with a wave vector K can be 
expanded in basis functions of a single irreducible rep- 
resentations of the space group of the crystal with wave 
vector K: 



e Ki/ 
B 0j 



(Al) 



where S^-" describes the spin vector of the jth magnetic 
ion in the 0th cell which determines the magnetic struc- 
ture of the crystal with a wave vector K, and So(^' ly |j) 
is the basis function transforming according to the z^th 
irreducible representation. The spin vector in the nth cell 
of the crystal can be derived from the spins in the zeroth 
cell by the equation 



(f U \3) 



„iK-t r . 



So(flj) 



(A2) 



In accordance with the Landau theory of second order 
phase transitions, the majority of magnetic structures 
are characterized by a single irreducible representation. 
However, there are cases that involve two or more ir- 
reducible representations of the magnetic space group. 
We will follow the method developed by Izyumov and 
NaishEj to calculate the basis functions, considering only 
magnetic structures defined through a single irreducible 
representation. 

The magnetic representation may be expanded into 
irreducible representations d K " of the wave vector group 
G K by 



(A3) 



"* = -77^ E Xm(9)x* K "(9) , (A4) 
U[Uk) hec° K 

xKg) = 6 h SpR h ^ ,K a ' ' M hu , (A5) 



where X Kl/ (<?) is the character of irreducible representa- 
tion d Kl/ of group Gk and \m is the character of the 
magnetic representation. Sh equals 1 for usual rotations 
and -1 for inversion rotations. SpR h is the trace of the 
matrix R^a of rotation for the group element g, which 



includes a rotation part h and a translation part Th- 



9*3 



= hi 



Th = Xj 



(A6) 



The summation in equation (A4) is taken only over the 
zero block of group Gk- The magnetic atomic compo- 
nents of the basis functions of the irreducible representa- 
tion for the 0th cell are 
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So(M*)= E S h dlfj((g)e-^^S, 



( < 



(A7) 
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TABLE III: Permutation of Cu atoms in La4CuaMoOi2 crys- 
tal by the elements of group P112i/m. a p is the returning 
translation vector. 



element 


atoms 
12 3 4 5 6 


cip 


{hi |000} 


12 3 4 5 6 


(000) 


{h 4 |00|} 


1 2 3 4 5 6 


(110)l-3, (111)4-6 


{h 25 |00|} 


4 5 6 12 3 


(110) 


{h 28 |000} 


1 2 3 4 5 6 


(OOO)i-a, (001)4-6 



TABLE V: Eigenvalues and eigenstates of the Hamiltonian 
H. The eigenstates listed in the table are not normalized. We 
use the following abbreviations: J = ^(Ji + J2 + J3), a = 
JiJ 2 -J!-J 3 A, b = Ji J3-JI-J2A, c = J 2 J 3 ~J?-JiA, A = 
Ji J2 - Jl + J3A, B = Ji J 3 - jf + J 2 A, C = J2J 3 - Jl + Ji A , 
and A = yj ' J\ + Jf + Jf - Ji J 2 - J 2 J 3 - J3J1 

E ^ 



(-3J-2A)/4 -2±£aP(i + ^/3a/3 + /3/3a 



TABLE IV: Irreducible representations of group C\ h for 
K=ia*. 



T4 




hi 


h 4 


h 2 s 


h 2 8 


Tl 




1 


1 


1 


1 


T2 




1 


1 


-1 


-1 


T3 




1 


-1 


1 


-1 


T4 




1 


-1 


-1 


1 


The 


space group of La4Cu3MoOi2 


is P112i/m. 


There 


are 


six 


magnetic 


copper ions in the 


primitive cell 


occu- 



pying the positions 

2(e) : l(af x , y x , 0), 2(x 2 , Vi, 0), 3(x 3 , y 3 , 0), 

4(1 -a*, 1 - 2/i 5 0.5), 5(1-012, 0.5), 
6(1 - ar 3 , 1 - 2/3, 0.5), 



+ £a/3a + f3aa 

1 (-3J + 2A)/4 -A&.app+cpap + ppa 

-^■aa/3 + %a/3a + /3aa 

2 3J/4 acta 

m 

a/3p + (3af3 + (3(3a 
aa/3 + a(3a + [3aa 

APPENDIX B: NEUTRON SCATTERING FROM 
A TRIMER 



where {x u Vi) are (1.0893, 0.8926), (0.8782, 0.4623), and 
(0.6465, 0.8816) for i=l, 2, and 3 respective^ The co- 
ordinates are written in the Kovalev systerrcil, and are 
shifted by (00^) from the coordinates listed in the Inter- 
national Table of Crystallography£3 

Group P112 i/m contains four elements, which are 
listed in table III. Tabic III also shows the permuta- 



tions of atoms by the action of the group elements^-Thc 
irreducible representations of the wave vector groupcil are 
listed in table ttVL The magnetic representation for the 
wave vector K=(^00) can be decomposed into irreducible 
representations as follows: 



l M 



3ri + 6r 2 + 6r 3 + 3t4 



(A8) 



The calculated basis functions are shown in table Q. 
Since the actual spin directions can be combinations of 
the basis functions within a single irreducible represen- 
tation, it follows from the table that the spins have to 
be either along c direction or in the a-b plane. For the 
irreducible representations T\ and t 3 , the spins on sites 
4-6 are antiparallel to those on sites 1-3, while they are 
parallel for the irreducible representations r 2 and T4. 



The Hamiltonian for a spin triangle with Heisenberg 
exchange interactions is 

H = JiSi • S 2 + J 2 S 2 • S 3 + J3S3 • Si . (Bl) 

The eigenstates and eigenvalues are listed in table |v|. The 
differential magnetic nautron cross section for inelastic 
transitions \S) -> \S') isB 



d 2 a 
dflcLu 



= Cqp(S) ( 5 a p - Q a Qp 



where 



Co 
P {S) 



a/3 



exp[iQ 



J2 (SM\Sf\S'M'){S'M'\S f j,\SM) 



MM' 



6 (huj + E(S) - E(S')) 



(B2) 



V I -^1) ^ F ^Q) exp [-2W(Q)} , 



m e c z 



Z _1 exp 



E(S) 



k B T 
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Also, Q = k — k' is the scattering wave vector, F(Q) is 
the magnetic form factorJH3 and Z is the partition func- 
tion. It can be shown that terms in the cross section with 
a =/= vanish. After averaging over all directions for Q 
we obtain the following cross section for a powder sample 



\ dflduj J . 



(4?-) 



\dfldu) ) x 



Cop(O)- [1+ 



sin(QR 12 ) (Ji 




-Ji) 


QR12 


A 2 




sin(QR 23 ) (Ji 


- J 2 )(J 2 


- J 3 ) 


QR-23 


A 2 




sin(QRi 3 ) (Ji 


- J 3 )(^3 


- J2) 



QR13 
5(huj + E - Ei) 



Cop(0)- 



sin(QR 12 ) 



QRi 



—a(a + b) 
i 2 + b 2 + ab 



sin(QRi 3 ) 
QR\3 
b(c - b) 



-be 



b 2 + c 2 - be 
ab 



b 2 



a 2 + b 2 + ab 
sin(QR 23 ) 



c 2 — be 
-b(a + b) 



QR23 

-c(c-b) 



6(huj - 
C p(l 



b 2 - 

- Eq 

2 



E2) , 
sin(QRi 2 ) 



QR 



12 



-A(A + B) 



-BC 



(B3) 



(B4) 



A 2 + B 2 + AB B 2 + C 2 - BC 



sin(QRi 3 ) 



AB 



QR13 \A 2 
B(C -B) 



B 2 + AB 
sin(QR 23 ) 



B 2 +C 2 - 


BC J 


QR23 


( -B{A - 






\A 2 + B 2 


+ AB 




-C{C- 


b) y 


X 


B 2 + C 2 - 


-BC) 



5(huj + E 1 -E 2 ) 



(B5) 



The spin triangles in La4Cu3MoOi2 are close to equilat- 
eral {R 12 = 3.39 A, i?i 3 = 3.50 A, R 23 = 3.41 A), so we 
can replace the Cu-Cu distances by R=i(i?i 2 + R 2 3 + 

Rl3)- 



d 2 a \ 
dildco I 



O-i-l 



sin(QR) 



QR 



S(Hu) + E a - Ex) 



d 2 a \ 
dtldcu J ^ 2 

d 2 a \ 
dildu J 



= Cap(0)-, 
5{huj + E 9 



1 



sin(QR) 



1 



QR 

E 2 ) , 
sin(QR) 



QR 



8{Tiw + Ei - E 2 



(B6) 

(B7) 
(B8) 



As each level is degenerate there is also an elastic cross 
section associated with intra-level transitions. However, 
this is not relevant for the experiment and so will not be 
listed. 
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